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Abstract. - Recent discovery of high T c superconductivity in Fe-based compounds may have 
opened a new pathway to the room temperature superconductivity. The new materials feature 
FeAs layers instead of the signature CuC>2 planes of much-studied cuprates. A model Hamiltonian 
describing FeAs layers is introduced, highlighting the crucial role of puckering of As atoms in 
promoting d-electron itinerancy and warding off large local-moment magnetism of Fe ions, the 
main enemy of superconductivity. Quantum many-particle effects in charge, spin and multiband 
channels are explored and a nesting-induced spin density-wave order is found in the parent com- 
pund. We argue that this largely itinerant antiferromagnetism and high T c itself are essentially 
tied to the multiband nature of the Fermi surface. 



Recently, a surprising new path to room-temperature 
superconductivity might have been discovered. The qua- 
ternary compound LaOFeP was already known to become 
superconducting below 7K [1], when its doped sibling 
LaOi-^F^FeAs (x > 0.1) turned out to have unexpectedly 
high T c of 26K [2]. Even higher T c 's were found by replac- 
ing La with other rare-earths (RE), reaching the current 
record of T c = 55K [3], These are the first non-cuprate 
superconductors exhibiting such high T c 's. 

The surprise here is that the most prominent charac- 
teristic of iron is its natural magnetism. By conventional 
wisdom, the high T c superconductivity in RE-OFeAs com- 
pounds is unexpected, all the more so since the super- 
conductivity apparently resides in FeAs layers. Follow- 
ing standard ionic accounting, rare-earths are 3 + , giving 
away three electrons, while As and O are 3~ and 2~, re- 
spectively. One then expects Fe to be in its 2 + configu- 
ration, two of its 4s electrons given away to fill As and O 
p-orbitals, with assistance from a single rare-earth atom. 
The remaining six d-electrons fill atomic orbitals of Fe in 
the overall tetragonal As/O environment of Fig. 1; the 
lower three t 2g orbitals should be filled while the upper 
two e g orbitals should be empty. However, the Coulomb 
interactions intervene via the Hund's rule: the total en- 
ergy can be reduced by making the spin part of the atomic 
wavefunction most symmetric and consequently the or- 
bital part of it as antisymmetric as possible, reducing 



thereby the cost of Coulomb repulsion. The simplest re- 
alization of this is to occupy a low t 2g orbital with one 
spin-up and one spin-down electron while storing the re- 
maining four electrons into the spin-up states. The result 
is a total spin S = 2 of Fe ++ , with the associated lo- 
cal magnetic moment and likely magnetism in the parent 
compounds. This is the situation similar to manganese, 
the Fe's nearest relative, whose five d-electrons feel the 
full brunt of the Hund's rule and typically line up into a 
large spin state, and very different from copper, where d- 
orbitals are either fully occupied or contain only a single 
d-hole, as in the parent state of cuprate superconductors. 
All told, the circumstances are hardly hospitable to any 
superconductivity, let alone a high temperature one. 

In this paper, we first argue that the above Hund's rule 
route to large local moment magnetism is derailed by sig- 
nificant banding effects, promoting enhanced itinerancy 
for most Fe d-electrons. We show how such itinerancy 
arises from the combination of two factors: a sizeable 
overlap among Fe and As atomic orbitals and the dis- 
tortion of the overall tetragonal structure into a locally 
near-tetrahedral environment for Fe ions, both generated 
by the crucial "puckering" of As atoms out of the FeAs 
planes (Fig. 1). The puckering rearranges the t 2g and 
e g crystal-field levels so that E t2 ~ E e - the situation 
"in between" the purely tetragonal (E t2g < E e ) and the 
purely tetrahedral (E t2g > E eg ) - thus bringing all d- 
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Fig. 1: (Colour on-line) (a) The three dimensional structure of 
the parent compound, with FeAs layer (Fe - red, As - green) on 
top of a REO layers (RE - yellow, O - blue); The blue square 
in the FeAs plane corresponds to the 'planar' unit cell (b). We 
denote two Fe atoms with A and B, while the two As atoms 
that are displaced up and down with respect to the layer are 
presented by dotted and crossed circles respectively. We give 
our choice of axes in the corner (note that some papers use a 
coordinate system rotated by 45 degrees), (c) The evolution 
of d-orbital energy levels from the tetragonal to tetrahedral 
crystal field environment. The puckering of FeAs planes results 
in the situation which is "in between", placing all d-orbitals 
near the Fermi level. 

orbitals into a close proximity of the Fermi level Ep, and 
maximizes direct overlap between Fe d- and As p-orbitals. 
The end result are numerous bands crossing Ep and a mul- 
tiply connected Fermi surface, containing both electron 
and hole sections. We introduce a two-dimensional tight- 
binding model which captures the relevant features of this 
multiband problem. Next, we argue that large number 
of broad bands and the absence of large local Fe moments 
betrays not only the failure of the atomic Hund's rule but, 
via the enhanced metallic screening, the absence of strong 
local correlations in general. This implies the key role for 
the nesting properties and we present an analytic calcula- 
tion of various responses for circular and elliptical bands 
relevant to this multiband problem. These responses allow 
us to account for the observed weak antiferromagnetic or- 
dering in parent materials and provide strong clues about 
the superconducting mechanism itself. In this sense, the 
Fe-based high T c superconductors differ from the hole- 
doped cuprates and are likely more closely related to either 
the electron-doped cuprates or other weakly to moderately 
correlated superconductors. 

The parent compound of the Fe-based superconductors 
has a ZrCuSiAs type structure [4], with eight atoms per 



unit cell, depicted in Fig. 1. The Fe atoms lie in a plane, 
same as O atoms precisely above them, in the adjacent 
REO layer. In contrast, the RE and As atoms (also lo- 
cated above each other) are puckered out of plane in a 
checkerboard fashion. The amount of puckering is sig- 
nificant, turning the in- plane tetragonal structure in the 
physically relevant FeAs layer into the nearly-tetrahedral 
one (the angle of the FeAs bond with respect to the ver- 
tical is 58.8° as compared to 54.7° for a tetrahedron [5]). 
As stated above, this has important consequences for pro- 
moting banding and rich orbital content near Ep. 

There available electronic structure calculations of 
LaOFeP [6], and of LaOFeAs, doped and undoped [7,8], 
consistently convey the key information that all five Fe 3d 
bands of are located at the Fermi level, in stark contrast 
with the cuprates. These bands are hybridized with 4p or- 
bitals of As/P located far below the Fermi level, centered 
between 6 and 2eV below Ep. There are five sections of 
the Fermi surface: two hole concentric, near-circular quasi- 
2D ones around the T point, two electron elliptical ones, 
centered around the M point, and a 3D hole band with 
a spherical Fermi surface around the Z point. Given the 
fact that the last one vanishes upon doping [8], and that 
the relevant physics appear to be two-dimensional, we will 
ignore this Fermi surface and neglect the interlayer cou- 
plings altogether. This idea is used in other works which 
aim to recreate the band-structure, either with all ten 
bands [9,10], or with some simpler minimal model [11,12]. 

To illustrate the key role of the puckering of As atoms 
on the electronic structure of FeAs, let us first consider 
the hypothetical situation in which all As atoms are pla- 
nar (Fig. 1). The tetragonal crystal field splitting pushes 
3d <2g orbitals (xy, xz, and yz) below the e g orbitals. 
In this arrangement, the overlap of Fe ti g orbitals with 
the neighbouring As p-orbitals either vanishes by sym- 
metry or is very small, the only source of broadening for 
these bands being the direct overlap of two d-orbitals on 
neighbouring Fe. The e g bands, on the other hand, do 
directly couple to the 4p-orbitals of As, but, since these 
p-orbitals are deep below the Fermi level, this coupling 
only pushes the e bands further up, increasing the crystal 
field gap. The consequence is that such a material should 
likely be a band insulator, turning into a local moment 
magnet once the Coulomb effects and the Hund's rule are 
turned on. A sizeable puckering changes the situation dra- 
matically: first, the Fe crystal field environment turns to 
near-tetrahedral instead of the tetragonal. In the purely 
tetrahedral case, the ti g orbitals (xz, yz, and x 2 — y 2 ) re- 
verse their position and are above the e g levels (xy and 
z 2 ). In the nearly-tetrahedral case of real FeAs layers, the 
tig orbitals are slightly above ei g , and the overlap due to 
the band formation makes all five bands important. This 
banding is the other crucial consequence of the puckering: 
since the p XtV orbitals are not entirely in the Fe plane, the 
overlap between these orbitals and xz, and yz d-orbitals 
becomes significant, and it actually provides the dominant 
contribution to electron/hole hopping. At the same time, 
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the hopping between x 2 — y 2 and p, or z 2 and p orbitals 
is only slightly reduced. 

Based on the above arguments, we construct a tight- 
binding model which incorporates the hoppings to the 
nearest neighbours and includes the relevant overlap in- 
tegrals. This model which reflects the key qualitative fea- 
tures of the electronic structure calculations [7, 8], and 
which can serve as the realistic platform for further ana- 
lytic calculations. As shown below, even such a simplified 
model must include all five d bands and is defined by the 
tight-binding Hamiltonian 

H = H + Ht+H int , (1) 

Ho = ^e<M a H a) +Y,e W P? )1 pf , ( 2 ) 

»»<* j,P 
it _ \^ . j(a)tJ/3) , 4-Fe j(«)t>') , 

where H describes local 3d and 4p orbitals, and H t ac- 
counts for Fe-As, Fe-Fe, and As- As hopping, in that order. 
Hi n t is the interaction term and will be discussed shortly. 
The operator d\ annihilates an electron in orbital a on 
Fe site i, and analogously, Pj annihilates an electron on 
site j in orbital (3. The summation over a takes into ac- 
count all five Fe 3d orbitals, but due to the doubling of 
the unit cell, there are actually ten of those bands, and 
the summation over j3 involves all three bands p x , p y , and 

Pz- 

The symmetry provides important guidance in under- 
standing the band structure of (1). For example, at Y 
point there are two doubly degenerate bands. One of these 
must be a symmetric combination - relative to A and B 
sites - of d xz , and d yz orbitals weakly hybridized with the 
As 4p bands, while the other is the antisymmetric combi- 
nation. The splitting between these two dublets originates 
both from the direct spread of the d xz / yz bands {t xzxz ), 
and from the p bands spread (tp s ). At the M point, these 
orbitals again two degenerate dublets, albeit in a different 
linear combination, which arc split only by the amount 
proportional to t F xz xz y From such analyses, we find the 
orbital energies and hoppings (all in eV's) 



a 


x 2 - y 2 


z 2 


xy 


xz 








£q 


-0.85 


-1.4 


-1.1 


-1.15 




X 


z 


t Fe 


-0.55 


-0.5 


-1.6 


-0.55 




-4.0 


-4.0 


^a.x/ y 


0.65 


-1.4 


1.5 


3.2 


j-As 
l 0,/3 


-0.8 


-0.45 


ta,z 


2.1 


1.25 




0.7 







The interband couplings are not tabulated and their 
values are t F f xy = 0.1, tl^ yz = -0.75, and = 0.8. 

The levels e a reflect our previous discussion of the crys- 
tal field splitting: E t . 2 ~ E e on the scale of t's. The 
hoppings reveal that the puckering of As atoms promoted 
d xz / yz bands to the physically most relevant ones, their 
coupling to the 4p orbitals being the strongest. These 
bands provide dominant content of the electronic states 




Fig. 2: (Colour on-line) The band structure (a), and the Fermi 
surfaces (b) following from the tight- binding Hamiltonian (1), 
and using the parameters of the tight binding fit. 

at Ep, where they get mixed with the other states, chiefly 
d xy , to finally form the two hole Fermi surfaces. Clearly, 
these effects are difficult to reproduce within a simple two- 
band model. In addition, significant mixing of different 
d-orbitals with opposite parity relative to the FeAs planes 
further reinforces the need to include the full d-orbital 
manifold into the basic description. 

Fig. 2 shows the band structure and the Fermi surface(s) 
following from (1). The key features of the Fermi sur- 
face are faithfully reproduced, with the central hole pock- 
ets nearly circular (actually, these are two ellipses which 
strongly interact and avoid crossing). In the vicinity of 
the M point, the two electron pockets have elliptical shape 
and do not interact at the crossing points located at the 
edges of the Brillouin zone. 

This brings us to Hi„ t in (1). The picture of puck- 
ered As atoms discussed above, promoting the bunching 
of local d-levels of Fe and their large overlap with As p- 
orbitals, indicates that the d-bands are near their optimum 
width, given the restrictions of dealing with 3d electrons 
and 4p levels far below of Ep. This reduces the impor- 
tance of the Hund's rule and points to the d-electron itin- 
erancy, rather than local atomic (ionic) correlations, as 
the most relevant feature. Indeed, this is consistent with 
the neutron scattering experiments [13], observing weak 
antiferromagnetism in the parent compound below 150K 
instead of the large local moment magnetism expected in 
the Hund's rule limit. The AF order is suppressed by dop- 
ing and ultimately vanishes in the superconducting state. 
This suggests that Hint should generically be comparable 
or smaller in magnitude than H t (1). For example, in the 
simple single band Hubbard model, with nearly circular 
(or spherical) Fermi surface, too large on-site repulsion U 
leads to the ferromagnetic Stoner instability, an itinerant 
prelude to the local moment formation dictated by the 
Hund's rule. We thus expect that the main effects of Hint 
can be understood by a detailed analysis of enhanced spin, 
charge and interband responses of the non-interacting part 
of H (1). 

With this aim in mind, we observe that various pock- 
ets of the Fermi surface depicted in Fig. 2 can be viewed 
as radial and elliptical distortions of the same idealized 
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circle, two of such ideal (hole) pockets located at T and 
two (electron) at M points. As long as such distortions 
are not too extreme, the responses in different channels 
can be evaluated analytically, thereby greatly facilitating 
theoretical analysis. Where comparison is possible, our 
analytic results appear to agree with numerical calcula- 
tions in Refs. [8,11,12]. 

We first look at the spin-susceptibility, and analyze how 
near-nesting of the Fermi surfaces can lead to SDW order. 
To do this, we have to make some assumptions about the 
Fermi surfaces and separate the most important contri- 
butions. While some nesting takes place within slightly 
deformed circular Fermi surfaces in Fig. 2, the main con- 
tribution to the enhanced response arises from similarly 
shaped hole and electron pockets, followed by a less im- 
portant one arising from different hole-hole and electron- 
electron pockets of the Fermi surface. This is easily ap- 
preciated by nothing that, for our idealized circles, the 
electron-hole nesting leads to a divergent contribution to 
the electron-hole propagator (i.e., an RPA bubble). So, 
we concentrate on the spin-susceptibility x s (q, lo) where 
one the particle propagators corresponds to the hole band 
at the r point with a circular Fermi surface and Fermi mo- 
mentum kp, and the electron band forming a slightly el- 
lipcitally deformed Fermi surface centered around M vec- 
tor, as depicted in Fig. 2. The electronic states at the 
Fermi level have kp(l + £) momentum if parallel to the M 
vector, and kp(l + rf) if perpendicular. The dispersions 
are 
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with m e ,h being the mass of the electron/hole band. For 
simplicity we assume that they have the same mass m e . 
Wavevector k is given relative to the M point in the 
case of the electron band, while 1 is defined with re- 
spect to the center of the Brilouin zone. Parameters 
£ and ij are tied to the eccentricity of the Fermi el- 
lipse as e — — 77 1 (2 + £ + 77)/ (1 + max(£,r))), and to 
the ratio of states enclosed by the two Fermi surfaces 
N e /Nh = (1 + £)(1 + 77). Below, we evaluate the particle- 
hole bubble due to the near nesting of only one hole and 
one electron band. Our results are universal, generally 
applicable to any situation involving elliptical Fermi sur- 
faces, and particularly relevant for FeAs, where one has to 
sum contributions due to nesting of each individual hole 
and electron band. 

If the eccentricity were zero, and the two bands had 
identical Fermi momenta (£ = 77 = 0), the real part of the 
susceptibility is would have simply been given by 



, . . m P A 

^ (q ,c = 0) = 2-log R 
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Fig. 3: (Colour on-line) The arrangement of Fermi surfaces 
with elliptical bands at the corners of the Brillouin Zone show 
in Fig. 2 (a), and the regularization of the singular part of the 
susceptibility due to the elliptical distortion of the electronic 
Fermi surface (b). For £ = 0, n — 0, the hole and the electron 
Fermi surfaces become identical and the susceptibility diverges. 

occurrs in Eq. (6) when q = M due to the perfect nesting 
of two hole and electron Fermi pockets. The nesting in 
FeAs is not perfect due to small distortions in Fig. 2, and 
this singularity is cut off. Still, it appears nevertheless 
that this particular response at q = M is the strongest 
incipient instability of our system. If Hi nt is overall repul- 
sive and not extremely weak, say modelled as a Hubbard 
repulsion on Fe sites, Udn 2 di , this instability will produce 
the spin density-wave (SDW) ground state at the commen- 
surate wavevector q = M. It seems natural to associate 
this Fermi surface instability with the observed weak AFM 
order of the parent compound [13]. 

To appreciate how the deformation of the electron Fermi 
surface cuts off the singularity, we now find the explicit 
expression for this more general situation. There are two 
different cases, depending on whether the two Fermi sur- 
faces intersect after one has been moved by M (so that 
their centers coincide). If they do not intersect (equiva- 
lent to £77 > 0), the susceptibility is set by 



x(,(q = M,a; = 0) =4 



m e (l +£)(! + v) 
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(7) 
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where A is the UV band cut-off. A logarithmic singularity 



where 3 = 1 + (1 + £) 2 and T = 1 + (1 + r/) 2 . 

Clearly, it is the last two term which cause the nesting 
divergence in the limit when the ellipse transforms to a 
circle (£ — > 0, 77 — > 0). When the Fermi surfaces do inter- 
sect (£77 < 0), the last two logarithms in Eq. (7) should be 
replaced by — log (2 + £ + ??) — log |5 — T|. This term is 
responsible for the singularity in this case. The divergent 
behavior of the real part of the susceptibility is shown in 
Fig. 3. 

Our analysis of the divergent part of the susceptibility 
was centered on the case when q = M, and the question 
remains whether that is the global maximum. The deriva- 
tives of the susceptibility with respect to q are well defined 
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due to the regularization by finite £ or v. It is trivial to 
demonstrate that the first derivative at q = M vanishes 
in all directions, which can alternatively be argued based 
on symmetry. Therefore, one has to look for the sign of 
the second derivative in both x and y direction in order 
to determine whether the susceptibility has a maximum, a 
minimum or a saddle point at q = M. Even if it turns out 
that the susceptibility has a maximum, it may be a local 
maximum, not the global one. While the general treat- 
ment of the problem will be presented elsewhere [14], we 
illustrate the situation by two circular Fermi surfaces with 
slightly different radii, kp, and fcp(l + £)- The susceptibil- 
ity due to the nesting of these Fermi surfaces is compared 
for the cases when q = M, and q = M + kp^n, with n 
being a unit vector pointing in an arbitrary direction. The 
former corresponds to concentric Fermi surfaces, the latter 
to two surfaces touching each other. The susceptibility in 
the former case follows as a special case of Eq. (7) 



xo(q = M,o) 
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and the result for touching circles is obtained by replacing 
the argument under the square root by 2. This is always 
slightly larger than the susceptibility following from Eq. 
(8), regardless the value of £. Such a result implies a dif- 
ferent ordering vector q = M + kp£_n, albeit only in a 
continuum theory. Our system is on a lattice, and M is 
commensurate with the reciprocal lattice, hence any insta- 
bility at that wave-vector will be enhanced by Umklapp 
processes, whereas this is not true for other incommensu- 
rate wave-vectors such as q. Furthermore, we may argue 
that two Fermi surfaces touching should not produce any 
unexpected divergences in the particle-hole channel, by 
simply observing the analytic expression Eq. (7) when £ 
or n = 0. 

Eq. (7) can be applied to all the possible pairs of 
hole/electron bands found in the band structure of FeAs. 
There are two circular hole surfaces of different radii paired 
with two electron surfaces which are the same, except that 
they are rotated by 90° - this just exchanges £ and r\ (the 
unmarked Fermi surface in Fig. 3a). For the UV cut off 
we choose the inverse lattice spacing. We now compare 
the relative values for the doped and parent systems, with 
the help from the band structure calculations. For the un- 
doped parent system, we estimate [6] £i w 0.27, r\\ w 0.45, 
£2 ~ 0.00, and 772 ~ 0.14, which yields x'o ~ 5.3m e at 
q = M. Doping moves Ep upwards, increasing the elec- 
tronic, and shrinking the hole Fermi surfaces. The corre- 
sponding surfaces are now further apart, so the susceptibil- 
ity is expected to be smaller. Using Ref. [8], we estimate 
£1 w 0.72, r?i w 1.11, £2 w 0.35, and n 2 w 0.65, giving 
Xo ~ 3.8m e at x = 0.1. Similar estimate is obtained by 
using our tight-binding band structure of Fig. 2. This is 
quite a bit smaller than the undoped value, and suggests 
rapid suppression of our SDW upon doping, as observed 
experimentally [13]. 



The SDW/AF order at q = M discussed above and 
observed in experiments, could in principle also be inter- 
preted in the local spin picture, as arising from the direct 
and indirect superexchange between Fe atoms. The direct 
superexchangc J\ is generated by the overlap of 3d or- 
bitals of neighbouring atoms, i.e., overlap between A and 
B atoms in Fig. lb. Two A(B) atoms, in contrast, have 
an insignificant direct overlap. However, from our band 
structure we know that bands d xz and d yz hybridize with 
4p orbitals of As. Let us for example take one A atom 
in a unit cell in Fig. 1, and consider its overlap with its 
next neighbour A to the right. Both of these atoms have 
their d xz orbitals hybridized with the Ap x orbital of the As 
atom standing between. The new hybridized bands both 
carry a significant fraction of the As orbital, so a hopping 
between these two atoms is enabled via the intermediate 
As atom. This hopping gives rise to the indirect superex- 
change coupling J2. Similar argument was presented in 
Ref. [15]. By the same mechanism, the indirect exchange 
between d x 2_ y 2 orbitals of neighboring iron atoms, due to 
their overlap with p orbitals of As, yields ferromagnetic 
nearest neighbor contribution to Ji [16,17]. Our earlier 
analysis suggests that the total Ji is significantly smaller 
than J 2 (C = |^i |/^2 -C 1). At such a high ratio of frus- 
trating AF couplings, the AF ordering takes place indi- 
vidually on A and B sublattices [18] irrespective of sign of 
J\. At the mean-field level, the relative order on the two 
sublattices does not affect the ground state energy since 
each B site interacts with four neighbouring A sites, two of 
these spins pointing in the direction opposite to the other 
two; consequently, there is no overall interaction. This im- 
plies that, on classical level, the ground state would have 
been degenerate with its ground state energy independent 
of the angle between two order parameters. Thus, we in- 
clude excitations - spin-waves - and investigate how their 
interaction affects the ground state. For this, we use the 
standard Holstein-Primakoff bosonization. Assuming that 
the angle between two order parameters on lattices A and 
B is 0, and introducing HP bosons a, and b on two respec- 
tive lattices, the following Hamiltonian is obtained 



H 



2SJ 2 ^{4(4a k + 4& k )- 

k 

(cos k x + cos fc y )(a k a^ k + & k ol k + h.c.) + 
k k 

2C(cos cos ^■)(aj c 6 k + 6[a k - alfi[_ k - a k &_ k ) 



2Ccos#(sin ^ sin ^)(a k 6 k + 6 k a k + a k &t k + a k fe_ k ) 

The Bogoliubov transformation of the Hamiltonian (9) 
yields two new excitations whose dispersions are given by 



ASJ 2 



(2 + cos k x + cos k y ± 2( cos -y- cos -^-) x 



k k 

(2 — cos k x — cos k y ± 2( cos 9 sin sin 



(10) 
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Fig. 4: (Colour on-line) (a) The zero point energy in units of 
JzSN as a function of angle 6 between two staggered magneti- 
zations, with £ = 0.5 [20]. The energy is the lowest for 8 = 0, n. 
(b) The staggered magnetization in units of ub per unit cell is 
plotted for three different spin values S = 1/2, 1, and 3/2 (from 
bottom to top, respectively). The dashed line corresponds to 
the experimentally observed value of 0.35/j.b [13]. 

We numerically determine zero point energy, and plot it 
in Fig. 4(a). The energy of the system is minimized when 
the spins on two sublattices are collinear in agreement with 
the experiments [19]. 

The staggered magnetization must be evaluated numer- 
ically for arbitrary f. Fig. 4(b) shows the magnetization 
(per iron site) due to quantum fluctuations as a function 
of £ for different values of spin S. The fluctuations are 
spin independent, but the resulting magnetization is not. 
Unrealistically low spin and large C are required in order 
to explain what is observed experimentally [13], thus hint- 
ing at significant itinerancy in the AF state, in line with 
our previous arguments. 

We now turn to the superconductivity itself, clearly the 
most difficult problem, ft is naturally tempting to use 
the above propensity for SDW at q = M in the parent 
compund to generate pairing interaction once the system 
is doped away from the AF order [21]. Following the ex- 
ample of electron-doped cuprates and various organic su- 
perconductors, this would imply an ordering of a nodeless 
kind, with electron and hole pockets in Fig. 2 fully gapped 
but with gap functions of different relative sign (see, for 
example, Ref. [8]). It is important to stress here the cru- 
cial role played by the Josephson-like interband scatter- 
ing between hole and electron Fermi surfaces in bringing 
about this form of superconductivity (see further below, 
and also [14]). Since hole and electron pockets are not 
identical, the gap magnitudes would not be either, but the 
difference could be quite small. Observing such a relative 
sign difference in the otherwise fully gapped state (an s± 
or s' state) would clearly be strong boost for this picture 
of superconductivity generated entirely by antiferromag- 
netic (SDW) spin fluctuations. However, the above change 
in sign implies sensitivity to ordinary (non-magnetic) im- 
purity disorder which could severely suppress T c and the 
gap. This effect, while still present on general grounds, 
appears less consequential in hole-doped cuprates, due to 
their strongly correlated, almost local nature. In Fe-based 
superconductors the correlations are not that strong, as 
we have just argued above, and this impurity sensitivity 
becomes an important issue. Finally, in order to generate 



the s± (or s') superconducting state, the interband inter- 
action - enhanced by the proximity to the SDW - must 
overcome the intraband repulsion (see below) , not an easy 
task [14]. 

In light of this, one should not out of hand dismiss 
the possibility that Fe-based superconductors are of en- 
tirely different kind from even the electron-doped cuprates 
and similar superconductors, where the purely repulsive 
interactions suffice to generate pairing near a magnetic 
instability. Their multiband nature could instead be a re- 
alization of the exciton-assisted superconductivity. The 
phonon interaction appears to be too weak to push T c to 
55 K by itself [8]. However, large number of highly po- 
larizable bands around the Fermi surface leads to strong 
metallic screening and a possibility of a dynamical over- 
screening, which turns fi, the familiar Coulomb pseudopo- 
tential of the Eliashberg theory, attractive in the certain 
finite wavevector and frequency regions. This would pave 
the way for the exciton-assisted superconductivity, long- 
anticipated but never unambiguously observed [22]. The 
basic idea is that the dynamical Coulomb interaction: 

v(q, U ) = ~r~~~^~\ , <q> w) = i + ^(q)xp(q, w) , (ii) 

(V(q) = 47re 2 /<7 2 ) turns attractive at some finite q and rel- 
atively low lo. Fe-based superconductors appear to have all 
the ingredients: their highly polarizable multiband Fermi 
surface produces neutral plasmon modes corresponding 
to electron and hole densities oscillating in phase (neu- 
tral plasmons). Such modes act as "phonons", particu- 
larly if m e and mh are sufficiently different. Furthermore, 
the nesting features lead to enhanced density response 
X P (q, w) near q = M and this could turn the effective 
interaction attractive at relatively low lo. Finally, the in- 
terband pairing [23] could still be essential, and acts to 
further boost T c irrespective of its sign: 

where X ee ihh) an d X e h are the e-e (h-h) and the interband 
coupling constants, respectively, and lo p is the character- 
istic frequency of the exciton-assisted pairing. T c gener- 
ated by this mechanism is notoriously difficult to estimate, 
both due to the competition from structural and covalent 
instabilities in the particle-hole channel and the need to 
consider local-field contributions to e(q, lo) [24]. Neverthe- 
less, this "hybrid" option - in which phonons help make 
intraband interactions attractive, or at least only weakly 
repulsive, and enable the magnetically-enhanced repulsive 
interband interactions provide the crucial boost to T c - 
should be kept in mind as the experimental and theoret- 
ical investigations of Fe-based superconductivity continue 
in earnest. 

In summary, Fe-based superconductors appear to offer 
a glimpse of a new road toward room-temperature super- 
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conductivity. We have constructed here a simplified tight- 
binding model which qualitatively describes the physics 
of FeAs layers where the superconductivity apparently re- 
sides. Analytical results were given for the elementary 
particle-hole response in charge, spin and multiband chan- 
nels and used to discuss various features of the SDW/AF 
order and superconductivity. We stress the importance 
of puckering of As atoms in promoting d-electron itiner- 
ancy and argue that high T c of Fe-based superconductors 
might be essentially tied to the multiband character of 
their Fermi surface, favorable to the s± (or s') supercon- 
ducting state. It is tempting to speculate that different 
T c 's obtained for different rare-earth substitutions might 
be related to the different amount of puckering in FeAs 
layers and regulating this amount might be the key to 
even higher T c . 

Acknowledgments. — The authors are grateful to I. 
Mazin, C. Broholm and C. L. Chicn for useful discussions. 
This work was supported in part by the NSF grant DMR- 
0531159 and by the DOE grant DE-FG02-08ER46544. 

Additional remark. — Since this manuscript 
was originally posted on the arXiv in April 2008 
(http://arxiv.org/abs/0804.4678) there were numerous 
significant experimental and theoretical developments in 
this exceptionally fast-paced field. Those most relevant 
to this work include the observation of the superconduct- 
ing gap in PCAR [25] and ARPES experiments [26-28] 
and various theoretical approaches exemplified by Refs. 
[29-35]. Especially notable among these are further theo- 
retical explorations of the s± (or s') superconducting state 
in Refs. [14,36-38]. 
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